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ABSTRACT 


This paper presents Reinforced Concrete (RC) slab design optimization 
technique for finding the best design parameters that satisfy the project 
requirements both in terms of strength and serviceability criteria while 
keeping the overall construction cost to a minimum. In this paper four 
different types of RC slab design named as simply supported slab, one end 
continuous slab, both end continuous slab and cantilever slab are optimized 
using three different metaheuristic optimization algorithms named as Genetic 
Algorithms (GA), Particle Swarm Optimization (PSO) and Gray Wolf 
Optimization (GWO). The slabs with various end conditions are formulated 
according to the ACI code. The formulated problem contains three 
optimization variables, the thickness of the slab, steel bar diameter, and bar 
spacing while objective involves the minimization of overall cost of the 
structure which includes the cost of concrete, cost of reinforcement and the 
constraints involves the design requirement and ACI codes limit. The proposed 
method is developed using MATLAB. Finally, to validate the performance of 
the proposed algorithm the results are compared with the previously 
proposed algorithms. The comparison of results shows that the proposed 
method provides a significant improvement over the previously proposed 
algorithms. 
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1. INTRODUCTION 

A reinforced concrete slab is a key structural component and 
is used in houses to provide flat surfaces (floors and 
ceilings). Concrete slabs are effective systems where putting 
columns interrupts the structure's (Audiences, parking lots, 
hotels, airports, etc.) serviceability to cover the lengthy 
spans. In recent past, metaheuristic optimization algorithms 
have been applied to many structural problems, and RC slabs 
are no exception as a result a number of articles on RC slabs 
optimization have been published. However, to find 
extensive use among practicing engineering, structural 
optimization algorithms need to be formulated as cost 
optimization and applied to realistic constructions that are 
subject to the real constraints of frequently used design 
codes such as the American Concrete Institute Code (ACI, 
1999). 


approach because the addressed optimization problem is 
highly nonlinear and multimodal and contains various 
complex constraints. In such cases sometimes, optimal 
solutions may not exist at ah. In most practical design 
problems, finding an optimal solution or even sub-optimal 
solutions is a difficult task. Furthermore, the metaheuristic 
approaches involve lesser computational complexity than 
deterministic approaches. 

The rest of the paper is arranged as the section II presents a 
brief literature review, in section III the mathematical 
modeling of different types of slabs are discussed, section IV 
provides a brief overview of metaheuristic algorithms used, 
section V presents the simulation results, and finally section 
VI presents the conclusions. 


In this paper the slabs with various end conditions are 
formulated according to the ACI code. The formulated 
problem contains three optimization variables, the thickness 
of the slab, steel bar diameter, and bar spacing while 
objective involves the minimization of overall cost of the 
structure which includes the cost of concrete, cost of 
reinforcement and the constraints involve the design 
requirement and ACI codes limit. 

In the presented work the metaheuristic optimization 
approach is utilized instead of deterministic optimization 


2. Literature Review 

K. C. Sharma et. al. [1], presented a review on cost 
optimization of concrete structures. These structures include 
beams, slabs, columns, frame structures, bridges, tanks of 
water, folded plates, shear walls, pipes and tensile 
components. It summarizes interesting and important 
outcomes and conclusions. The paper also includes a review 
of cost optimization based on reliability. The results of such 
research will be of great value to engineering practitioners. 
S.M.R. Tabatabai et. al. [2], introduced a system for optimal 
reinforcement design and non-linear reinforced concrete 
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structure analysis. The specially tailored ORCHID (Optimum 
Reinforced Concrete Highly Interactive Dimensioning) 
program is used for the design and optimization of 
reinforcement. F. Ahmadkhanlou et. al. [3], presents a 
general formula for cost optimization of single-and multi¬ 
span RC slabs with different end conditions (simply 
supported, one end continuous, both end continuous and 
cantilever) subject to all ACI code constraints. M.G. Sahaba et. 
al. [4], presents the cost of optimizing reinforced concrete 
flat slab buildings under the British Code of Practice 
(BS8110). The objective function is the building's total cost 
including floor, column and foundation costs. The cost of 
each structural element is that of reinforcement, concrete 
and shaping material and labor. A genetic algorithm is used 
for a global search in this hybrid algorithm, followed by a 
discrete form of the method Hook and Jeeves. B.A. Nedushan 
et. al. [5], This article deals with cost optimization of single¬ 
way concrete slabs in accordance with the latest American 
Code of Practice (ACI 318-M08). The goal is to minimize the 
slab's total cost including concrete and reinforcement bar 
costs while meeting all design requirements. Particle Swarm 
Optimization (PSO) is used to solve the restricted problem of 
optimization. As PSO is designed for unconstrained problems 
of optimization, a multi-stage dynamic penalty was also 
implemented to solve the constrained problem of 
optimization. A. Kaveh et. al. [6], The cost optimization of a 
single-way reinforced concrete floor system consisting of a 
hollow slab is presented in this article. The system cost is 
considered the objective function and the design is based on 
the ACI 318-05 standard of the American Concrete Institute. 
This function is minimized using the harmony search 
algorithm, subject to design constraints. A. Kaveh et. al. [7], a 
new modified particle swarm optimization algorithm (PSO) 
is used in this paper to optimize the design of large-scale 
pre-stressed concrete slabs. The modification is 
accomplished by adding some probabilistic coefficients to 
particle velocity and is called probabilistic swarm 
optimization of particles (PPSO). These coefficients provide 
the algorithm with simultaneous exploration and 
exploitation and decrease PSO's dependence on its 
constants. A. H. Gandomi et. al. [8], a new metaheuristic 
optimization algorithm called cuckoo search (CS) for solving 
structural optimization tasks is introduced in this study. 
Combined with Le'vy flights, the new CS algorithm is first 
verified using a nonlinear benchmark restricted optimization 
issue. CS is subsequently applied to 13 design problems 
reported in the specialized literature to validate against 
structural engineering optimization problems. The CS 
algorithm's performance is further compared with different 
algorithms that represent the state of the art in the area. For 
the most part, the optimal solutions obtained by CS are far 
better than the best solutions obtained through existing 
methods. A. Akin et. al. [9], This paper presents the 
application of the harmony-based search algorithm to the 
optimal detailed design of special seismic moment 
reinforced concrete (RC) frames under earthquake loads 
based on American Standard specifications. The objective 
function is selected as the total frame cost that includes the 
concrete, formwork and reinforcing steel costs for individual 
frame members. M. Aldwaik et. al. [10], a model for cost 
optimization of reinforced concrete (RC) flat slabs of 
arbitrary configuration in irregular high-rise construction 
structures is presented in this article. The model is general 
and can include any combination with or without openings 
and perimeter beams of columns and shear walls in the 
plane. For flat slabs of arbitrary configurations, a general 


cost function is formulated taking into consideration not 
only the cost of concrete and steel materials but also the cost 
of construction. Using Adeli and Park's robust neural 
dynamics model, the nonlinear cost optimization problem is 
solved. The methodology has been applied in a real-life 36- 
story building structure to two flat slab examples. Not only 
does the methodology automate the RC slab design process, 
it also results in cost savings of 6.7-9 %. 

3. Model Formulation 

This section presents the mathematical modeling of slab 
designs and formulation of the objective functions. 

3.1. One-way reinforced concrete slabs 

3.1.1. Objective function 

A total cost function can be written as follows: 

C t = C c +C r + C f} #( 1 ) 

Where C a C* and Cf are, respectively, the costs of concrete, 

reinforcing bars and formwork and finishing materials. For 
any given location, the cost of formwork does not vary 
significantly and can therefore be dropped from formulation 
(Ahmadkhanlou and Adeli [3]). The definition is as follows: 

C c = LbhCm 2 ) 

C r = w 3 LA s (^#(3) 

Where L, b, h, Cj-, A s> and C * are the span length, span 

width, slab thickness (Fig. 1), concrete cost per unit volume, 
unit weight per unit volume of steel, cross section area of 
reinforcement bars, and reinforcement bar cost per unit 
weight, respectively. The calculation of quantity A 3 is as 

follows: 


ndz fb\ 


Where and s are the reinforcement bar diameter and 
spacing, respectively. 
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Figure 1: Typical cross-section of RC slab. 

3.1.2. Design constraints 

As mentioned earlier, cost function optimization is based on 
the limitations defined by ACI 1999 code [29]. The 
constraints included flexural constraint, shear constraint, 
limitation of serviceability, and limitation of deflection. As 
stated below, they are defined and expressed in a normalized 
form. 

FLEXURAL CONSTRAINT 

Nominal flexural strength, 'M n , should be greater than the 
ultimate design moment, 

M 

*W = nr-is»,* = t9, #(S) 

c|>M n 
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In Eq. (4.5), is calculated as follows: 

^=^,#(6) 

Where and k are the clear span length and the moment 
coefficient for the continuous slab depending on the slab 
support type, respectively. Table 1 shows the values of k. In 
Eq. (4.6) the maximum moment coefficient value given in 
Table 2 for four different support conditions (simply 
supported, continuous at one end and simply supported at 
the other end, continuous at both ends and cantilever). In Eq. 
(4.6), w is a uniformly distributed charge factored. The 
AhmadkhanlouandAdeli article [8] considered loading cases 
as suggested by the 1999 code of the ACI [30]: 

w = 1.4 x (DL x & x DL S } + 17 x LL x fc, #(7) 

Where DL, LL, and DL 3 are dead floor loads excluding slab 

self-weight, live load, and slab self-weight. The calculation of 
DL 3 as follows: 

DL 3 - (bh - A 3 )w c + A 3 w sr #(8) 
where W c is the weight of the concrete per unit volume. 
Table 1: Moment coefficient for continuous slabs 


Exterior Span 

Support 

Middle 

Support 

i 

~24 

i 

+ I4 

i 

“15 

Interior Span 

Support 

Middle 

Support 

i 

“IT 

i 

+ 16 

i 

“u 


Table 2: Maximum moment coefficient, k, used for design 
of RC slabs. 


Simply 

One End 

Both End 

Cantilever 

Supported 

Continuous 

Continuous 

i 

i 

i 

i 

8 

1C 

IT 

2 


The nominal bending moment, Af„, is calculated as follows: 
Where f y is the specified yield strength of the reinforcement 

bars and the corresponding depth of the concrete 
compressive stress block from which it is calculated (Fig. 1). 
AJy 


= 


k- m 


Q.S5 fjb 

where f c is the specified compressive strength of concrete. 

SHEAR CONSTRAINT 

The nominal concrete shear strength, should be greater 
than the ultimate factored shear strength, V u ; 

*aOO = ! £ D, * = Q.8E, #(10) 


The ultimate factored shear force is defined as follows: 

wl n „ v 

Vu = K-f, #(n) 


where k^ is the shear coefficient for continuous slab that 
depends on the type of slab supports. The values of k v are 

given in Table 3. The nominal shear strength of concrete is 
defined as follows: 

V c = 2 J&'bd, #(12) 

Table 3: Shear coefficient for continuous slabs. 


Simply 

Supported 

One end 
Continuous 

Both ends 
Continuous 

Cantilever 

i 

1.15 

i 

2 


SERVICEABILITY CONSTRAINT 

The percentage in one-way RC slabs of longitudinal 
reinforcement steel, p, and bar spacing, S, should be between 

minimum and maximum limits permitted by the design 
specification. 


0 a CO = 


■-IS 0, #(13) 


Pmax 

0400 = —— ISO, #(14) 

p 

0 5 M = ^-1S&, #(15) 

0&CO = —- 1 ^ Oj #(16) 

3 max 

where the p max is given by: 

Pmax = Q-75pb> #(17) 

p b is defined as follows: 

„ U { 87000 \ , v 

fc - MSft f(e75o+4j'* (l ® 

fc{ 800 \ 

-^zhccTZ) 


in S, L units,p b ■ 


where f c f and f y are in MF&#(18a) 
and is calculated from 

forfl < 4000psi J ^ : = Q.85, 

"ft = 

f 0.SE for f! ^ 4QQQpst 


, ffj - 4QQQ\ 

I 0.85 - 0.05 100Q J > 0.65 for fc > 4000ps£ 


###(19) 


in S.I.units,for fj S 27.S6 MPa, fa = 0.BE, 

ft = 

0.05 for fc < 27.50Aff > 3 

(0.85 - 0.05 ( /l? ~^ 3 ) > 0.65 forfc > 27.58 


The minimum area of flexural reinforcement longitudinal) is 
selected as follows: 


[0 QOZO&ii, for steel graded C and 50 
10 0 018 & ft, for steel grad es SO 


#(2fl) 


and the minimum and maximum bar spacing are defined as 
follows: 

5™, = max(r,d & ),#(21) 

in S. Limits, S min = max (25.4^1^ d t ) J #(21a 
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S max = min (18", 3ft), #(22^ 

inS.Lunits, = min (45 /.2mm, 3/i), #(22a) 

DEFLECTION CONSTRAINTS 

Slab thickness, ft, shall not be less than the minimum slab 
thickness, ft 

= ^-1<C,#(23) 

Where fthas a minimum thickness of 1.5 in (38.1 mm) in 

Table 4. Table 4 values apply to normal concrete weight and 
fy = 413.7 MPa (60,00Q psi)- For f y other than 

413.7 MPa (60,000 psi), % (specified in Eq. (24)) 
multiplies the values. For lightweight concrete with w € 
between 14.14 feitf/m 3 (SOI.fr//t a ) and 

18.06 kN/m? (llSi.fe//f 3 ), the values must be multiplied 
by (as set out in Eq. (25)). 


4. Column offset in the direction of offset from either axis 
between centerlines of successive columns shall not 
exceed 10 percent of the span. 

5. All loads shall be distributed uniformly over a whole 
panel due to gravity only. 6. Unfactored live load shall 
not exceed the unfactored dead load twice. 

6. Unfactored live load shall not exceed the unfactored 
dead load twice. 

These constraints include flexural restrictions, one-way 
shear constraints, two-way shear constraints, serviceability 
constraints, and deflection constraints. The limitations are 
explained below and expressed in a standardized form. 

FLEXURAL CONSTRAINT 

The design moments are distributed across each panel in the 
direct design method. As shown in Fig 2, the panels are 
divided into columns and middle strips. In each strip, 
positive and negative moments are achieved. The column 
strip is a slab with a width on each side of the centerline of 
the column. 


a ± = 0.4 + ■ 


fy 


#(24) 


fy 


100,000 

in S.Lunits, a ± = 0.4 -I- ^, 

where f y in AfPa#(24a) 

= ma^l.bS - O.OOSw,, 1.09), #(25) 

Table4: Minimum thickness for solid one-way slab 
according to ACI code 


Simply 

One end 

Both ends 

Cantilever 

Supported 

continuous 

continuous 

L 

L 

i 

L 

20 

1A 

28 

To 


3.1.3. Design variables 

Single-way RC concrete slab design variables consist of three 
variables: slab thickness (ft), reinforcement bar diameter 
(ii b ), and reinforcement bar spacing (s). Slab thickness and 

reinforcement spacing can be considered as integer 
variables, such as multiples of 5mm and 10mm respectively 
in the SI system, or multiples of 1/tf " or 1/4" in the usual US 

system. Since it is necessary to assign the diameter of the 
reinforcement bars from limited numbers, it must be 
considered as a discrete variable. ACI supplies eleven 
different bar sizes with a diameter of Q.375"(9.53 mm) 
from bar size #3 to bar size #16 with a diameter of 2.257" 
(57.33 mm). 


b c = min(^, ^),#(2 6) 

& m = ! 2 -2x& e .,#£27) 


where h c and h m are the width of half column strip and 
middle strip, respectively. 

coEihim iUsp 



Figure 2: Middle and Columns Strips. 

Positive nominal bending strength at the middle and 
negative nominal bending strength at both ends of the 
column strip and middle strip, should be greater than 


3.2. Reinforced concrete flat slabs 
3.2.1. Design constraints 

This paper analyzes the concrete flat slabs by Direct Design 
Method based on [29]. For using direct design method, ACI 
presents six limitations; 

1. In each direction, there have to be at least three 
continuous spans. 

2. Successive span lengths (center-to-center supports) 
must not differ by more than one-third of the longer 
span in each direction. 

3. The panels must be rectangular, with the ratio of panel 
dimensions longer to shorter, the center-to-center of 
supports measured, and not exceed 2. 


the ultimate design moment in column strips (column) and 
middle strips (mid) at those ends, 

9iM = ^z~1 <D,#(2B) 

M + 

ISO, #(29) 

SscM = - 1 S 0, #(30) 

M" 

SimM =^Z~l<0, #(31) 

M + 

SsmM = - 1 < a #(32) 
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expressed in Table 5. and Table 6. M c is the total static factor 
for a span defined as: 


ffSnM = ^ - 1 < o, #(33) 


In Eqs. (4.28-4.33), = 0,9, and Af u is calculated as follows: 

M ucs = 34) 

= (1 - (35) 

Where k m is the distribution of the total span moment 

coefficient for four different span conditions (inner span, 
outer edge unrestrained, outer edge fully restrained, slab 
without beam between inner support and edge beam) and k 
is the portion of the negative inner moment, the negative 
outer moment, and the positive moment resisted by the 
column strip. The values of k m and k a are respectively 


M,= 


8 


#(36) 


Where In is the clear range of supports in the direction of 
considered moments, the face-to-face measurement of the 
supports is not less than G.GSJi, and the factored uniformly 
distributed load is defined as: 

w u — 1.2 x (DL + + 1,6 x LL, #(37) 

Where is the reinforced concrete density and h the slab 
thickness. 


TableS: Distribution of total span moment, k m . 



Exterior edge 
Unrestrained 

Slab without beam between 
interior support and edge beam 

Exterior edge 
fully restrained 

Interior 

span 

Interior negative 
factored moment 

0.75 

0.70 

0.65 

0.65 

Positive factored 

moment 

0.63 

0.52 

0.35 

0.35 

Exterior negative 
factored moment 

0 

0.26 

0.65 

0.65 


Table 6: Portion of interior negative moment, exterior negative moment, and positive moment resisted by column strip, k Q 



h 

h 


0.5 

1 

2 

Interior negative moment 

0.75 

0.75 

0.75 

Exterior negative moment 

1 

1 

1 

Positive moment 

0.6 

0.6 

0.6 


The nominal bending moment, is defined as follows: 

M ncc ,l = 4 C3 fy (d - |), #(38) 


M. 


HFTli & 


= 35) 


Where a is the corresponding depth of the concrete 
compressive stress block calculated as follows: 

^scotfy 


^cal — 


0.85// x 2h c 
A, 


#{40} 


_ ‘amidfy . . 

“ 0,85 fibj #(41) 


the quantity of ,4, is calculated by; 

A Sfol = '^i^+ l),#(42) 


In effective slab width, a fraction of factored slab moment 
resisted by the column, should be less than the 

nominal flexural strength, 

g 7 (x) = iso,<j> = os, #(+4) 

The effective slab width shall be the width of column plus 
l.Sh of slab. M sc in interior column (int col) and edge 
column (edge col) and yy are defined as follows: 

= omquu + O.Sfl JlalJ- q' D X(m W 

where q' D , l' 2 , and l ^ refer to the shorter span 
edge col ~ 0.3M C ,#(46) 
i 


Yf =' 




,#(47) 
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The maximum values for yy is provided in Table 7. 

Table 7: Maximum modified values of yy . 


Column location 

Span direction 


[within 

Maximum modified y f 

Corner column 

Either direction 


2 : n,no4 

1.0 

Perpendicular to the edge 


2:0.004 

1.0 

Edge Column 

Parallel to edge 


2:0.010 

1.25 

-r* 1 ,n 

-— 1.0 
»<djt 

Interior column 

Either direction 


>0.010 

1.25 

- = < 1.0 

‘♦(Dj; - 


Where b ^ is the shear perimeter length perpendicular to the 
bending axis, and is the shear perimeter length parallel to 
the bending axis, c± is also the width of the column 
perpendicular to the bending axis, whereas c 2 is the width of 
the column parallel to the bending axis. Fig. 3 calculates 
these perimeters. 

U * ** .1 <t Col p 


In Eq. (4.52), is defined as follows: 


D| c IA 


c ? + d 


r 


C3 k. 


v u.CD t 


Critical 

-- T —--j section 


—f 



u.AB 


Shear 

stress 


W interior column 


Ci + c#2 


c 2 4 


4 Col 


-f-r 







AB 


a? 

Shear 

stress 


Critical 

_i j section 

C C D I ^ 

C 

ft) Edge column 

Figure 3: Assumed distribution of shear. 


ONE-WAY SHEAR CONSTRAINT 

The nominal concrete shear strength, should be greater 
than the ultimate factored shear strength, V u ; 

03 OO = - 1 =S Q, <J> = 0.75, #(48) 

The ultimate factored shear force is calculated as follows in 
the interior spans (int) and edge spans (edge): 

W = ^i 2 (y-rf),#(49) 

/l.l Sl„ \ 

Vu B ^ B = ^uh[-^- d J^( 5 °) 

The nominal concrete shear strength is defined as: 

^ -K- 

TWO-WAY SHEAR CONSTRAINT 

The shear stress strength, should be greater than the 
ultimate factored shear stress, V u ; 

g 9 (x) = -^-l<04 = 0.75, #(52) 

Sk V n 


^=^ + ^,f ( S3) 

A c Jc 


Where V U£ = the ultimate shear calculated by Eq. (4.55), 
A c = concrete area calculated by Eqs (4.55-4.56) along the 
assumed critical section. f c = the assumed critical section 
property analogous to the Eqs (4.57-4.58) calculated polar 
moment of inertia, c = the distance in the critical section 
calculated by Eqs (4.59-4.60) between the central axis and 
outlines, and in Eq (4.61), fv is given. (4.61): 


V ug = w fl (l t l 2 -b 1 b 2 '),#0>*) 

A cinteol = 2{hi + b 2 ) x d, #(55) 


Ar.arfjar.nf — (2 £>i + ij) X d, #{56) 

fbl bid 3 . v 

hinteol = d [f + ^ 1 J + ^S7) 

= d ^ - (2 bi - b 2 ) x c z j + #(S8) 


/, 


c edge col ' 


C int col = y.'#(59) 

b l 


^oigoool ~ 2^+bz 

Yv = i -rrMGty 


. #(60) 


Where C± is the column width perpendicular to the bending 
axis and is the column width parallel to the bending axis. 
For two-way members without shear reinforcement, the 
shear stress strength is calculated by: 

v n =v c = min ^4^ - 2^, #(62) 


In Eq. (62), for interior columns, the value of a 3 is 40 and for 
edge columns is 30. 

SERVICEABILITY CONSTRAINS 

The area of reinforcement bars, should be greater than 
the minimum area of reinforcement, -4 am i n , and bar 
spacing, S should be between minimum and maximum limits 
allowed by the design specification in reinforced one-way 
slabs. 


S 106 O = : 


-- 1 < 0, #4,63 
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S 11 OO =— -ISO, #4,64 

5 

tfia(r) = — - 1 < a 

Table 8 presents the minimum area of flexural 
reinforcement: 


Table 8: The minimum area of flexural reinforcement, 

4 


Reinforcemen 
t Type 

fyPSl 


Deformed 

bars 

<413.7 MPa 

0.002 

Deformed 
bars or 
welded wire 
reinforcement 

>413.7* (Pa 
[somopsi) 

f 60000 \ 

max ^0.0013 X ——— A M ,€JQQ14A 3 J 

iTtS.Ltsmts r f x in MPa 

f 413.7 \ 

max j 0.001S X -y—A sr 0.00144 j J 


The minimum and maximum spacing of the bar is defined as 
follows: 

*min - ",d bl -d esg j, #{G6] 

in S.L units, s min = 

d b , | , #( 66 s) 

where d a g is the diameter of aggregate. 

min (18 \ 2 ft), at cri ti cal secti ons 

min (1B% 3ft), at other sections * ^ ^ 

inS.I.units, s m „ x = 
min (4S7 2mm, 2ft), at critical sections 

m in (4b / 2mm, 3 ft), at o ther secti ons ^ a * 

in this paper the formula for critical sections is assumed for 
all sections of the RC flat slabs. 

DEFLECTION CONSTRAINS 

Slab thickness, ft, shall not be less than the minimum slab 
thickness, ft 

*lbOO = ^-l<Q,¥(68) 

where is presented in Table 9. 


Table 9: Minimum thickness of slabs without interior 
beams, h min . 


fyrPSi 

Without drop panels 

Exterior Panels 

Interior Panels 


Without edge beams 


40000 

In 

in 

33 

36 

60000 

in 

in 

3fi 

33 

75000 

in 

in 

28 

31 


where l r , is the clear span in the long direction, measured 
face-to-face of supports. 


3.2.2. Design variables 

The concrete compressive strength, the slab thickness 
(ft), the reinforcement bar diameter (dj, and the 
reinforcement (5) spacing were included as the design 
variables. The diameter numbers of the reinforcement bars 
(d b ) vary for four end spans because they differ in the 

ultimate design moment. The thickness of the slab (ft) and 
the spacing of the reinforcement (s) may be considered as 
integer variables (multiple of 5mm and 10mm respectively), 
while discrete variables must be assigned to the compressive 
strength of the concrete [fj) and the diameter of the 

reinforcement bars (d b ). 


Table 10:List of possible values for fj. 


S. No. 

Concrete strength (£,') 

1. 

14MPfl (2000psi) 

2. 

t7MPa (2 EOOps) 

3. 

21 MPa (3000ps0 

4. 

24 MPa (3500 psi} 

5. 

2S MPa (40C0psi) 

6. 

31 MPa (4EC0psi) 

7. 

34MP (5000 psi) 

8. 

41 MPa (6000psi) 

9. 

5E MPa (SOCOpsi) 

10. 

69 MPa (lOOOOpsi) 


Table 11: List of possible values for . 


S. No. 

Diameter of reinforcement bars («vL ) 

1. 

10mm 

2. 

12mm 

3. 

16mm 

4. 

20mm 

5. 

25mm 

6. 

28mm 

7. 

32mm 

8. 

36mm 

9. 

40mm 

10. 

50mm 


3.2.3. Objective Function 

A total cost function can be defined as follows for the 
reinforced concrete flat slab: 

C t =C^C c +C r +C f lH69) 

Where C c is the concrete cost. C c can be computed as: 

C ff = ^^#(70) 

Where ij. is the span length, center to center of supports in 
the direction in which moments are considered, i 2 is the 
span length, center to center of supports in the transverse 
direction to and ^ is the cost ratio of the unit volume of 
concrete to the unit volume of concrete C f and C) 

are the costs, respectively, of negative and positive 
reinforcement bars in interior and exterior supports, 
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shaping and finishing materials. As mentioned in section 3, it 
is possible to drop the cost of formwork from the 
formulation. 

4. OPTIMIZATION ALGORITHMS 

The design goal in optimizing a design could be simply to 
minimize cost of production or maximize efficiency in 
production. An optimization algorithm is a process that is 
performed iteratively by comparing various solutions until 
an optimal or satisfactory solution is found. With the advent 
of computers, optimization has become part of computer- 
aided design activities. The metaheuristic optimization 
approach is used in the presented work instead of 
deterministic optimization approach because the 
optimization problem addressed is highly nonlinear and 
multimodal and contains different complex constraints. 
Sometimes there may be no optimal solutions at all in such 
cases. Finding an optimal solution or even sub-optimal 
solutions is a difficult task in most practical design issues. In 
addition, metaheuristic approaches involve less complexity 
computational than deterministic approaches. 

4.1. Genetic Algorithm (GA) 

GAs uses probabilistic transition rules to give research a 
way, as opposed to numerous methods. Creating a tool that 
depends on mutually exclusive events such as throwing a 
coin is neither a simple arbitrary hunt nor a choice. 
However, GAs use random choice as a tool to accompany a 
search to the search space areas with expected progress. In 
the following stages, the basic principles of a GA are 
explained. 

INITIAL POPULATION 

Population of randomly selected chromosome-shaped 
solutions. 

TERMINATION CRITERIA 

In practice, this is the termination criterion; for instance, if 
the current generation t > t max the maximum permissible 

number of generations is given, then the end condition is 
given. 

OBJECTIVE OR FITNESS FUNCTION 

IT defines the goal the GA needs to achieve. 

SELECTION 

Genetic algorithm reproduction begins with selection. 
Selection is a method of selecting separate strings 
(chromosomes) for next generations. 


ClLTOilLUSilOK 

I Gi L i 

- -- ■ ■— 



Figure 4: Presentation of Gene, Chromosome and 
Population in GA. 

CROSSOVER 

Crossover is nothing but a reproductive process that creates 
different individuals in successive generations by combining 
bits from two chromosomes of the previous generation. 


MUTATION 

Mutation is the process that haphazardly alters genetic 
information to avoid trapping into local minima. 

NEW GENERATION 

We get new chromosomes at the end of the generation 
through production processes such as selection; crossover 
and mutation applied to a population of n chromosomes 
until a new set of n persons is created. Thus, this set 
becomes the new population. 



Figure 5: Flowchart of Genetic Algorithm (GA). 


Particle Swarm Optimization (PSO) 

PSO is a population-based stochastic optimization 
technology based on the social behaviors observed in 
animals and insects, e.g. bird flocculation, fish schooling and 
animal herding [27]. 

BASIC STRUCTURE 

In PSO, through a systematic approach, a number of particles 
are moved in search space. Each particle i have a position, 
JCj(t), and a velocity, (f), in PSO, at time t. The particles 

current position and best position ever are stored in a 
memory. The particle velocity will be changed based on the 
historical data stored in the memory as well as random 
information. The new velocities will be used to update the 
particles' current position, 

INITIAL POPULATION OR SOLUTIONS 

A number of initial solutions are needed for the proposed 
algorithm to initiate space exploration solution. Basically, 
these initial solutions are the particles used in the search. 
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PARTICLE POSITION, VELOCITY, UPDATION AND STOPPING 
CRITERION 

Because L particles are used by the proposed PSO algorithm 
to explore the feasible space. Initially, the PSO algorithm 
generates L random numbers as particle velocity to update 
the particle position. After updating iteration in for each 
particle, the appropriate velocity ensures the feasibility of 
each particle. The algorithm also needs to modify particle 
velocities during the search to guide the particles through 
the region's most desirable regions. Originally, the PSO 
algorithm uses the following equation to update the speed: 
^k = Vk + ClTlf-PiitifljC- — ^k) + ^2*2 — -Pk). # 

Where V k is the particle at k th velocity, is the current 
particle (solution) at k m iteration, Pfetejt * s the particles 
personal best, * s the particles global best, TJ and ^ is 

a random number between (0 and 1) Assume 
ri = n.7R,r7 = 0.4R. C u C 2 are learning factors (or) social 
and cognitive parameters. Usually equals to C 2 and 
ranges from 0 to 4. Assume C-_ = C 2 = 1. 

With the new velocity \ t the particle's updated position 

equation is given as: 

F k -P k +V k r# 

The two above equations imply that a new design is 
searched for the global optimum by using the velocity vector 
being explored based on the local and global bests. Thus, 
through neighborhood learning and previous design, the PSO 
design update is done. 



Figure 6: Flowchart of PSO Algorithm. 


GREY WOLF OPTIMIZATION (GWO) 

The Grey Wolf Optimization (GWO) was proposed by 
Mirjalili et al [28]. The GWO is inspired by the social 
structure and hunting behavior of the gray wolves. 

The gray wolves are strictly following the social hierarchy of 
leadership. The alpha (a) wolf, still at the top of the 
hierarchy, leads the hierarchy group. Likewise, after alpha, 
beta (P) wolf is called the wolfs second level, the wolves of 
the third and fourth level are called delta ((3) and omega ( to) 
respectively. The alpha wolf is followed by all wolves (beta, 
delta and omega), while delta and omega follow the beta 
wolves, and delta wolves follow only omega. Since the omega 
remains at the lowest level, they have no followers. 

Hunting is now led by alpha, beta and delta wolves and 
restful wolves (omega). The movement of the entire 
population is guided by the top three best solutions in the 
optimization problem and these solutions are referred to as 
alpha, beta and delta respectively the other solutions are 
considered as omega. 

ENCIRCLING THE PREY 

The first hunting step is to surround the prey. Gray wolve's 
encircling process is equivalent to the optimal solution being 
encircled by all the population and is given by: 

D = \C-X pr „ y (i) - X waif td\,#Z.B 
■WO + 1) = Xpreyit) ~ T D, *3.9 

- -i 

Here i represent the current iteration number, while A and C 
are the coefficient vectors, X woi f and X prvy are the wolf 

- -i 

and prey position vectors respectively. The A and C vector 
coefficients are calculated as follows: 

A = 2a-r 1 -a, # 3.10 

C = 2r 2 #3.11 

In eq. 3.10 the value of vector a is linearly decreased from 2 
to 0 with the iterations and r 2 are random vectors within 
the interval of [0,1]. 


tM 



Figure7: The position updating process in GWO as 
presented by Mirjalili et al [4]. 

HUNTING 

The prey position is known in the actual hunting scenario, 
but the optimal solution is not known in the problem of 
optimization, so alpha, beta and delta solutions estimate a 
rough estimation of the optimum location knowing they have 
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the best knowledge of the solution. The update of the wolve's 

position is as follows: 

-#■ -#■ — 


-j- 

Xyfj tf 


o +1) = 


x 1 px 2 +x 3 

3 


i# 


The above equation estimates how the wolf (population) 
should move to obtain the prey (optimum) and the mean of 

-r -r 

the possible locations of pray. The X u X% and X% are the 
expected position of prey (optimum solution) on the basis of 
the position of ff, and S respectively. These positions are 
estimated as follows: 

X l = X a - A ± ■ (D a l. D a = \C 1 -X a - x waif \ r # 

X 2 = Xg — A 2 • (i3p}, Dp = | Ci ■ X@ — X wo if |j # 
X E =X S - Az ■ (D s ), D S = \C 1 -X S -X waif l# 


ATTACKING 

As the gray wolf starts tightening grip to prey the movement 
of the prey becomes smaller and smaller as the wolves move, 
and finally the prey stops moving and the wolf performs the 
final attack. In a mathematical model, the vector value (a) to 
be reduced from two to zero by every iteration in the linear 
movement of prey and wolf (population position) is reduced 
and finally, the prey (optimum) gets the movement of the 



Figure 8: Flowchart for Graywolf Optimization 
Algorithm (GWO). 


5. Simulation Results 

Optimization findings are provided in this chapter for four 
kinds of one-way RC slabs with distinct support 
circumstances. Table 12 presents the prevalent data used for 
the simulation. The initial design variables values are 
selected according to the ACI code for the design of concrete 
slabs. is set at $'13QQ/£cm (StASfkg) for 

reinforcement steel. The concrete price differs with the 
strength of the concrete as shown in Table 13. Practical 
values for variables h and s are assumed to be E mm and 
10 mm multiples respectively. 


Tablel2: Common Data used for Simulation 


Parameter 

Parameter 

Parameter 

Name 

Description 

Value 

fy 

Specified yield 
strength of 
reinforcement bars 

275,0 MPa, 
or (40 ksi ) 


Weight of steel per 
unit volume 

76.372S kNfm 3 
or (4yui&//-t a ) 

fc 

Specified 

compressive strength 
of concrete 

ZU.6B MFa 
or (3 ksi) 


Weight of concrete 
per unit volume 

23.6kN/m 3 , 

Cover 

Cover 

19,05 mm, 
or(0,75 in) 

c\ 

Cost of concrete per 
unit volume 

99,5 

nr(7fiS/nyri.y 

L 

Span length 

3,96 m, 
or (13 ft) 

DL 

Dead load of floor 
excluding the self¬ 
weight of slab 

0,48 KN/m 2 , 
or (10 lb/ft 1 ) 

LL 

Live load 

2,39 KN/m 2 , 
or(i0lbfft 2 ) 

Cl 

Cost of 

reinforcement bars 
per unit weight 

1300 S /ton, 
or (1,43 $/kg) 


Tablel3: Concrete specified compressive strength fr 
and its cost, cK 


14 MPa (2000psi) 

_“C_ 

93,655 S/m 3 

17MPa (2 5 OOps) 

96,94 S/m 3 

21MPu (300Ops 0 

99,56$/m 3 

24 MPa (3500psi) 

102,18 S/m 3 

2BMPa (4000psi) 

106,765 S/m 3 

31 MPa (4500psi) 

108,73 S/m 3 

34 MP (EOOOpsi) 

110,695 S/m 3 

41MPfl (SOOOpsi) 

126,415 S/m 3 

55 MPa (SOOOpsi) 

206,98 S/m 3 

69 MPa (lOOOOpsi) 

293,44 S/m 3 
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Table 14: Genetic Algorithm (GA) Configuration. 


Parameter's Name 

Parameter's Value 

Population Size 

100 

Max Generations 

200 

Initialization Technique 

Random (Uniformly 
Distributed) 

Crossover Technique 

Scattered 

Crossover Probability 

0.8 

Mutation Technique 

Random (Gaussian 
with mean = 0) 

Mutation Probability 

0.01 

Selection 

Stochastic uniform 


Table 15: Particle Swarm Optimization (PSO) 


Configuration. 


Parameter's Name 

Parameter's Value 

Number of Particles 

100 

Max Iterations 

200 

Initialization Technique 

Random (Uniformly 
Distributed) 

Inertia Weights 
([^‘m in ■” ^Wi-a ]) 

[0.1000 1.1000] 

Acceleration Coefficients 

([ci^aD 

[1.49,1.49] 


Tab le 16: Gray Wolf Optimization (GWO) Configurat ion. 


Parameter's Name 

Parameter's Value 

Number of Wolfs 

100 

Max Iterations 

200 

Initialization 

Random (Uniformly 

Technique 

Distributed) 


ism 

iSOO 

1400 

10M 

3 

^ iom 

■J 0 

| ** 
600 

400 


47 2502 Henn: 47.K0-J 


Sop i 3 s*-5t ES^npr^ticHi 

Figure 9: Convergence Plot of Genetic Algorithm (GA). 

Pflflt Fundkm VaW; 47 z&n 


lli-v-nlmH 


Figure 10: Convergence Plot of Particle Swarm 
Optimization (PSO). 


BgeE f itnns 5 = 4-J.^EDU 



Figure 11: Convergence Plot of Gray Wolf Optimization 
(GWO). 


Table 17: Optimization results for simply supported slab. 



Slab thickness 


Rebar Spacing (:.) 

Total Cost (!;!j! 

Ref [31 

171.45mm [6.75111] 

9.525mm (0.375in) 

165.1mm (6.5in) 

26.45 

Ref [51 

158.75mm (6.25in) 

12.7mm (0.5in) 

228.6mm (9.0in) 

26.57 

Ref [23J 

158.75mm (6.25in) 

15.875 mm (0.625in) 

368.3mm (14.5in) 

26.36 

GA 

160.02mm (6.30in) 

12mm (0.47in) 

350.52mm (13.8in) 

23.61 

PSO 

160.02mm (6.30in) 

12mm (0.47in) 

350.52mm (13.8in) 

23.61 

GWO 

160.02mm (6.30in) 

12mm (0.47in) 

350.52mm (13.8in) 

23.61 

Table 18: Optimization results for one end continuous slab. 

uMlldsuliE 

Rebar Spacing (:.) 

Total Cost (Sb) 

Ref [31 

141.48mm [5.57in] 

9.525mm (0.375in) 

177.8mm [7.0in] 

22.98 

Ref [51 

133.35mm (5.25in) 

9.525mm (0.375in) 

139.7mm (5.5in) 

22.76 

Ref [231 

133.35mm (5.25in) 

12.7mm (0.5in) 

254.0mm (lO.Oin) 

22.78 

GA 

135.13mm (5.32in) 

12mm (0.47in) 

398.78mm [15.7in] 

20.06 

PSO 

135.13mm (5.32in) 

10mm [0.39in] 

289.56mm (11.4in) 

19.89 

GWO 

135.13mm (5.32in) 

10mm [0.39in] 

289.56mm (11.4in) 

19.89 
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Table 19: Optimization results for both ends continuous slab. 



Slab thickness 

Rebar diameter 

Rebar Spacing ( : ) 

Total Cost (\) 

Ref [3] 

120.65mm (4.75in) 

9.525mm (0.375in) 

177.8mm (7.0in) 

19.93 

Ref [5] 

114.3mm (4.5in) 

9.525mm (0.375in) 

139.7mm (5.5in) 

20.64 

Ref [23] 

114.3mm (4.5in) 

12.7mm (0.5in) 

254.0mm (lO.Oin) 

20.5 

GA 

115.06mm (4.53in) 

12mm (0.47in) 

289.56mm (11.4in) 

19.1 

PSO 

115.06mm (4.53in) 

10mm (0.39in) 

340.36mm (13.4in) 

16.95 

GWO 

115.06mm (4.53in) 

10mm C0.39in) 

340.36mm (13.4in) 

16.95 


Table 20: Optimization results for cantilever slab. 


Slab thickness Rebar diameter Rebar Spacing (:.) Total Cost 


Ref [3] 

342.9mm (13.Sin) 

9.525mm (0.375in) 

50.8mm (2.0in) 

60.22 

Ref [5] 

317.5mm (12.Sin) 

15.875mm (0.625in) 

317.5mm (12.5in) 

59.31 

Ref [23] 

317.5mm (12.Sin) 

17.145mm (0.875in] 

241.3mm (9.5in) 

59.96 

GA 

320.04mm (12.6in) 

10mm (0.39in) 

119.89mm (4.72in) 

47.33 

PSO 

320.04mm (12.6in) 

16mm (0.63in) 

309.88mm (12.2in) 

47.25 

GWO 

320.04mm (12.6in) 

16mm (0.63in) 

309.88mm (12.2in) 

47.25 



Rrsnlraxomfin^ bar? taping in mrtars 

Figure 12: Variation of the total cost function versus s 
while the other two design variables h and d are kept 
constant. 



Figure 14: Variation of the total cost function versus 
d_b while the other two design variables s and h are 
kept constant. 





n.l 015 US r>25 0.3 095 0.$ CMS 

fhlcKnes* gl (h > in 

Figure 13: Variation of the total cost function versus h 
while the other two design variables s and rf L , are kept 
constant. 


Figure 15: Variation of total cost function versus h and 
s when d b is kept constant to 0.0160m. 


- s = D 21 jnrf = DJdUi 
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IKK 

1LKL 



Figure 16: Variation of total cost function versus . & 
and s when h is kept constant to 0.3200m. 



Figure 17: Variation of total cost function versus 
and h when s is kept constant to 0.3100m. 


The cost optimization model described in the paper can also 
be used as a tool for conducting parametric research, gaining 
insight into the entire design, and gaining trend data and 
answering questions when asked. For instance, Fig. 12-14 
displays variation in the overall cost function versus a single 
design variable h, ii b , or 5 while the other two design factors 
remain fixed. Fig. 12 depicts that the cost reduces very 
quickly when bars spacing (£) is increased up to 100mm, 
however the further increase in bar spacing doesn't affect 
the cost that much. On the other hand, The Fig. 13 shows that 
cost increases linearly with rise in slab thickness (ft). Finally, 

the Fig. 14 shows that the cost rises drastically with rise in 
the size of the reinforcement bar (i 

Fig. 15-17 displays variation in the overall cost function 
versus two design variables while the third design variable 
kept fixed. Fig. 15 depicts that the complete cost function 
variation versus h and s if is held constant. This figure 

indicates a significant cost increase when the spacing of the 
bar is low (less than 80mm). Fig. 16 Displays the complete 
cost function vs. and s if h is held constant. This figure 

indicates a significant rise in cost function with reduction in 
the size of the reinforcement bar and increment in spacing of 


bar. Fig. 17 Displays the complete cost function vs. h and 
if s is held constant. This figure indicates a significant rise in 
cost function with rise in the size of the reinforcement bar 
and rise in slab thickness. 

6. Conclusions 

In this research, cost optimization of RC slabs with distinct 
support conditions (simply supported, one end continuous, 
both end continuous and cantilever) was provided using the 
three distinct optimization methods (Genetic Algorithm 
(GA), Particle Swarm Optimization (PS0), and Gray Wolf 
Optimization (GWO)). The slab design was based on ACI 
code as discussed in section 3, and the procedure involved 
finding the optimum slab thickness, reinforcement bar 
diameter, and reinforcement spacing the details of these 
variables are presented in section 3. 

According to the outcomes, the GWO, which was first used to 
optimize the concrete slab, showed (refer to Fig. 9-11) better 
convergence speed to optimize concrete structures relative 
to the GA and PSO algorithms. The GWO achieves the best 
result in just 15 iterations (Fig. 9), while the PSO and GA take 
around 25 iterations (Fig. 10 and Fig. 11). Furthermore, the 
GWO and PSO both achieve similar and better outcomes than 
the GA for one end continuous slab (Table 18Error! 
Reference source not found.), both end continuous slab 
(Table 19) and cantilever slab (Table 20), although for 
simply supported slab all three algorithms give similar 
results. 

The comparison of the presented method with the 
previously proposed methods, shows (refer to Table 17-20) 
that the presented method achieves better outcomes than 
the previously proposed works. The presented method 
achieves 10.43% cost reduction (Table 17) for the simply 
supported condition, 12.61% cost reduction (Table 18) in 
one end continuous support condition, 14.95% cost 
reduction (Table 19) in both end continuous support 
condition and 20.33% cost reduction in cantilever support 
condition (Table 20). 
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